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ABSTRACT 

Stauffer (1973) developed a theoretical model for the prediction 
of air flow above ocean waves. In this study his results were extended 
with the addition of buoyance effects and two different lower boundary 
conditions. 

Through nximerical solutions it was possible to determine that: 

1) mean wind velocities at lower levels closely approximate the initial 
mean velocity profile for large wave numbers (k ^ 0.45m or small 
wave amplitudes, 2) mean wind velocities at lower levels fluctuate more 
from the initial mean velocity profile as the wave number decreased or 
the wave amplitude increases, 3) solutions were very sensitive to the 
level where the velocity gradient is computed in the lower boundary 
condition, 4) disturbance potential temperature at lower levels in- 
creases inversely with wave number, and 5) the presence of temperature 
fluctuations has an insignificant effect on the Reynolds stress or on 
the stream function. The numerical solutions show general agreement 
with the observational analyses of wave modified wind profiles by 
Davidson (1974) . However, these comparisons are not conclusive because 
of the sensitivity in the lower boundary condition. 
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I 



INTRODUCTION 



The forerunner of this study (Stauffer, 1973) used a 
theoretical numerical model to investigate wave-related 
momentum transfer. The model was used to obtain the Rey- 
nolds wave stress of the perturbation-veloci ty field in 
the atmospheric layer above ocean waves. The study was 
similar to that of Yefimov (1970) and utilized an initial 
value numerical procedure similar to that of Newman (1969). 
Equations were solved by utilizing central differencing in 
general but with the turbulence terms evaluated at the pre- 
vious time step. Solution separation was avoided by in- 
corporating the Matsuno (1966) finite difference scheme 
every 50 time steps. Experiments were run varying the in- 
dependent variables: wave number (k) , dynamic velocity 

(U.), roughness parameter (z ), and the coefficient of 
turbulent exchange (V). The measure of the error in the 
change of the stress in air (R) was computed at each time 
step and calculations continued until the stress reached 
a steady state condition defined by R ^ 1 x 10 Maximum 

Reynolds stress occurred at about one meter above the wave 
surface, and the largest momentum transfer was found to 

exist for the largest mean wind speeds, lowest wave numbers, 

2 

and a turbulent coefficient of approximately 0.02m /sec. 

The magnitude of the Reynolds stress increased as the 
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mean wind velocity approached that of the wave phase 
speed in the lowest levels of the boundary layer. Thus 
when the wave velocity is close to the mean wind velocity 
in the lower levels^ the disturbance energy is the 
greatest . 

The emphasis in this study is on the changes in the 
mean wind which are caused by the’ wave stresses which 
were computed by Stauffer (1973). His basic model is 
modified to include buoyancy effects which arise during 
non-neutral conditions. Also the lower boundary condi- 
tion was changed to include the mean wind shear. Two 
approximations to this boundary condition were investi- 
gated. Experiments were carried out in which the follow- 
ing quantities were varied: wave number (k) , grid size, 

convergence criteria, wave amplitude (a) , and lower 
boundary condition. It was found that the mean wind de- 
parted from the logarithmic profile in the lower region 
when the wave number became small and when the wave ampli- 
tude was large. Some of these variations were consistent 
with the observations of Davidson (1974). 
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II 



DEVELOPMENT OF THE MODEL 



Development of the model is similar to that of 
Stauffer (1973). The Boussinesq equations for two dimen- 
sional motion take the following form; 



3t 



+ W . VW = - Vir - 
3 



g 0 k 

0 



+ |f , 



■|| + W • V0 = 



3u 3w 
3x 3z “ ° ' 






( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 



where 



° * h (-u'W) , 



'■z = ^ ^ > ' 

1C 

IT = 0 (P^/P) + gz r 

p o o 

0 = T (p /p)^ - 0 , and 

o o 

K = R/c 

P 

1C 

Here T (p /p) is the total potential temperature, 0 is 
o o 

the mean potential temperature for the whole region, R 
is the universal gas constant, is the specific heat 

at constant pressure, x is the horizontal axis, and z is 
the vertical axis. Equation (2.1) is the equation of 
motion; equation (2.2) is the first law of thermodynamics, 
and equation (2.3) is the continuity equation. The 
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quantities u*, w*, and 0* are turbulent fluctuations. A 
straight line above a variable signifies the stationary 
part (time averaged) of each variable, and a wavy line 
represents the non-s tationary part (space averaged) . The 
velocity (W) and the potential temperature ( 0) are 
broken up as follows: 

W = (U(x,t) + u(x,z,t)) |i + w Ik , and (2.4) 

0 = 0(z,t) + 6(x,z,t) . (2.5) 

Following Yefimov (1970 ) , the turbulent fluxes are approxi- 
mated as follows; 



u’w' 




-w ' 6 ' 



-u'6 ' 



8z 



T = K + V ( 



, 8 u 

8z 



8 w 

ax' 



|2 

w' , 

^ ae ^ . ae 

^0 az ^0 az ' 



V, 



ax * 



( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 



Here K and are the coefficients of turbulent exchange 
u 

for the mean fields, and V and are the coefficients of 

u 

the wave-caused fluctuations. 

The following vorticity equation is formed by taking 
|j • Vx with equation ( 2 . 1 ): 



a? ^ „ a? ^ a^u a? ^ ac 

3 — + U -r— + w :r + u X— + w x— 

. dt dx N 2 dx 3z 



a^ , au ,au 9w,, , a^u 

_ [K + v(x— + X— ) ] - V( + 

^2 dz dz dx '\2_ 

3z 3x 3z 



^3 
3 w 

3x^ 



)- 



g_ ^ 

0 3x ' 
o 



(2.10) 
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where 



r - -Lii 

8z 3x 



(2,11) 



In this derivation, equations (2.3), (2.4), (2.5), (2.6), 

and (2.7) have been used. When equations (2.4), (2,5), 

(2.8), and (2.9) are introduced into equation (2.2) it 



becomes 



39 ^ ,39 ^ 30 ^ 39 ^ 39 

3t ^3x '^3z “3x '^3z 



~ 2 2 
9_(k ^)+v (i-i + i-® 



) 



(2.12) 



These derivations have used the following assumptions; 



3t 




0 



(2.13) 



Take the x-average of equations (2.10) and (2.12); 



3? 3^ 

^3x '^3z 






30 30 

^3x '^3z 



3 , 30, 

^ ^^0 



(2.14) 



(2.15) 



Subtract (2.14) and (2.15) from (2,10) and (2.12), respec- 
tively, and neglecting products of the fluctuation quanti- 
ties gives : 



ii 

3t 



+ 




w- 



2 

3 U 



2 3 3 

3 , , 3u . 3w, , , 3 u . 3 w, 

..2 dz dx -N 2. „ 3 

3z dx dz dx 



g ^ 

0 3x ' 

o 



li 

3t 



+ 



,.30 







3^0 

3x^ 



3^0. 

2^ 

3z 



(2.16) 

(2.17) 
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since the motion is non-divergent (equation (2.3)), a 
stream function can be defined such that: 



u 





(2.18) 



When these relations are introduced into (2.16) and (2.17), 
they become: 



3 ^ ^ ^ _ 321 






3x‘ 

tv(- 



- + + li i!u 

3x3z^ 3x^ 3z^ 



3z 9z 9x 9x 9z 



d ip g ^ 



9x 



0 9x 
o 



li = + M 10 + u ri!i + illi 

3t ~^3x 3x 3z 2 _ 2^ * 

9x 9z 



(2.19) 



( 2 . 20 ) 



In this study the motions are forced by a surface wave of 
wave number (k) and phase speed (c). Therefore, the stream 
function and the potential temperature departure are 
written as : 

if; = A(z,t) cos [k(x-ct)] + B(z,t) sin[k(x-ct)] ,(2.21) 

6 = D(z,t) cos [k(x-ct)] + E(z,t) sin[k(x-ct)] .(2.22) 



Substitute relations (2.21) and (2.22) into (2.19) and 
separate the cosine and sine terms. Equating the coeffi- 
cients of the cosine terms gives: 



|^[-^ - k^A] = k(U-c)(Bk^- 1-^) + Bk 1-| 



3z 

+ — „[v(^ 

3z 3z 



+ Ak ) ] + Vk [ 



3z 
2 ,3^A 



3z 



„ 2 ' 0 
3z o 



(2.23) 
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-M 












The coefficients of the sine terms give: 






9z^ 

2 



- k B] = -k (U-c) (Ak - 



9^A 



)- Ak 



9z‘ 



9z 9z 



+ Bk ) ] + Vk [ 



2 , 9^6 



+ Bk ] + 



9^U 

9z^ 

gkP 



9z‘ 



0 



(2.24) 



If the same procedure is carried out with equation (2.20), 
the cosine coefficients give: 



9d 

9t 



= -k (U-c) E + 



1 ® 

9z 



kB - 



VgtDk 



9z^ 



(2.25) 



The sine terms give: 



If = -k(c-U)D + II kA - VQ[Ek‘ 



l!l] 

9z^ 



Substitute (2,21) into (2,14) yields 






9^B 

9z^ 



- B 



9^A 

9z^ 



] = 



9z 






(2.26) 



(2.27) 



Equation (2.15) will not be used in this study because it 
turns out that the temperature effects are very small. 
Therefore 0 is equal to its initial value in this study. 
The equations (2.23) through (2.27) form a complete set 
in the variables D, E, and U. 
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III. BOUNDARY CONDITIONS AND INITIAL CONDITIONS 



Yefimov (1970) defined the boundary conditions for the 
lower boundary of the turbulent layer as: 



U(Tl,t) + = -akc cos [k(x-ct)] , 

w (x,U/t) = akc sin [k(x-ct)] , 

where the height of the wave surface is given by: 



(3-1) 

(3-2) 



0 = a cos [k(x-ct)] . 



(3.3) 



1/2 



The deep water wave phase speed is: c = (g/k)^'“ . (3.4) 

If (3.1) and (3.2) are in terms of the stream function 



(2.18), they can be written as: 

U(rift) + (x,r)ft) = akc cos [k(x-ct)] , 

djp 



(3.5) 



- ■ 3 ^ (x,ri/t) = akc sin [k(x-ct)] . 



(3.6) 

The stream function and its derivatives at z = r] can be 
approximated by the values at z = 0 , but U(T]ft) cannot be 
approximated by U(0,t) since U has a very rapid change 
near the surface. Therefore U(U/t) is expressed with a 
Taylor series expansion as: 

(3.7) 

Substitute (2.21) into (3.5) and use (3.7) and (3.3) which 



9U 

u(n,t) = u(o,t) + ^ n • 



gives : 



If (0,t) = -a[kc . If] , 



8b 

dz 



( 0 ,t) = 0 



(3.8) 

(3.9) 
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The condition U(0,t) Z 0 has been used. When (2.21) is 
substituted into (3,6), the following are obtained: 

A(0,t) = ac , (3,10) 

B(0,t) = 0 o (3,11) 

The following boundary condition is used for the potential 
temperature fluctuation : 

0(x,O,t) = 0 (3,12) 

This condition is not as accurate as the velocity conditions 
above, but the results shown later indicate that the tem- 
perature effects will be small in any case. Substituting 
(2,22) into (3,12) yields: 

D(0,t) = E(0,t) = 0 (3,13) 



In order to close the problem, the following conditions are 
imposed at the top of the region: 



A(H,t)= -||^(H,t)= B(H,t)= -||(H,t)= D(H,t)=E (H,t)= 0 ,(3.14) 



where z = H is the upper boundary. 

The initial conditions are given by: 

A(z,0) = ac sinh [k ) ] /sinh kH^ 

B (z , 0) = 0 , 

D (z , 0) = 0 , 

E(z,0) = 0 , 

and U(z,0) = ^ » 



, (3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 
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1 




where 



roughness 



z = 
o 



: = H-Az, = dynamic velocity, 

parameter, K = 0.35 = von-Karmen*s constant, and a = wave 
amplitude. The equations give an irrotational initial 
state with no temperature fluctuations. The initial mean 
wind is logarithmic. 

The following fields which are independent of time 
take the following form: 



K = K U^z , 

- T^ 

e = i„ (f_ * 1) , 

o 



(3.20) 

(3.21) 



where T^ is the temperature scale. 
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IV 



NUMERICAL PROCEDURES 



Centered finite differences are introduced into equa- 
tions (2.23) through (2.26) with the exception that the 
turbulence terms are evaluated at the previous time step. 

The left hand side of the equations are solved for the 
time tendencies with the Guass elimination technique which 
is described in Richtmyer (1957). Centered time differenc- 
ing is used except that the solution is restarted every 50 
time steps with the finite difference scheme developed by 
Matsuno (1966) . This procedure eliminates solution 
separation. 

Note that the value of U given in equation (3.19) with 
the coefficient given by (3.20) leads to a zero value for 
the right hand side of (2.27). Therefore U can be expressed 
as follows : 



U(z,t) 




1 ) 



+ u (z , t ) . 



(4.1) 



Substitute (4.1) into (2.27) and 

with respect to z which yields: 

z 



u 




dz 



9 



integrate three times 



(4.2) 



o 

• 1. ,9b, 3a„, , . , , 

where T = -uw = - — k (-^A - — B) (4.3) 

2 dz d z 

is the Reynolds wave stress. In this derivation the right 
hand side of (2.27) was rewritten and the coefficient K was 
replaced by the constant V . 
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The boundary conditions (3,8), (3,9), (3.13), and 



(3.14) are applied in a straightforward manner with a 
one-sided derivative when required. 

The time integration of the system of equations is 
accomplished by first solving the finite difference form 
of (2.23) through (2.26), and then computing the new U 
from equation (4.1). Integration is continued until a 
steady state is reached. The following quantity is used 
as a measure of the steadiness of the solution; 



The integration is continued until R > R# where R will 

m — m 

take on three different values. 




1 

2 



R 





(4.4) 



i = l 
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V 



NUMERICAL SOLUTIONS 



The following constants were used in all numerical in- 
tegrations : 





U* 


= 0 




Z 


= 7 




o 






V 


= 0 


Three 


forms of 


the 


these 


numerical 


int 


lowing 


condition 


wh 



0.167 m sec ^ , 

X 10 m , 

. 2 -1 



(5.1) 



l£ 

dz 






-akc 



(5o2) 



The 
singula 
tion, P 
can be 
level z 



direct use of (3.19 
rity when applied at 
rofessor R. Davis ha 
alleviated by comput 
= z^. In this case 

3A 
3z 



) in equation (3.8) 
z = 0. In a priva 
s suggested that th 
ing the derivative 
the boundary condi 



U. 



( 0 , t ) = -a [kc + 



K(z + z ) 
1 o 



leads to a 
te communica 
is problem 
at a false 
tion becomes 

(5.3) 



This formula reduces to (5.2) when z^ is very large. This 
procedure does not take into account changes in the mean 
wind near the surface through u. 

A third boundary condition was used which changes 
as the mean wind changes. This takes the form: 
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4 ' 

<4 



... afSi 

m 



I 






« 





3a 

32 



(0,t) = -a[kc + 



U(Az,t) - U(0,t) 
Az 



] (5.4) 



Numerical solutions were obtained varying six indepen- 
dent parameters as shown in Table I. These were wave 
number (k ) , wave amplitude (a ) , measure of the error in 

the change of the stress in air (R ). false level (z,), tem- 

m 1 

perature scale (T ) , and the boundary condition -^(0,t). 
Experiments were carried out as exhibited in Table II for 
the different values shown in Table I. 

Figures 1 and 2 show the distribution of mean wind 
velocity (U) with height (z) for series 1 and 2, respective- 
lyo These experiments use the lower boundary condition 
used by Stauffer (1973). Both figures show that the mean 
wind is close to the logarithmic profile for the largest 
wave numbers. They also show large variations as the wave 
number decreases. These departures are principally in the 
lower layers below approximately six meters. The smallest 
wave number (k = 0.050 m ^) shows a much larger velocity 
than the logarithmic profile with a first maximum at approxi- 
mately 1.5 meters. The character of the solutions for the 

two values of R are similar, but for particular wave 
m 

numbers there are considerable differences. 

Figures 3 and 4 represent series 3 for wave numbers 
k =.0.35 m ^ and k = 0.20 m respectively. These actually 
represent the Reynolds wave stress (T) as a function of 
height. Curve six in each figure is for the case z^ = 
infinity which corresponds to the boundary condition used 
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TABLE !• Values of the variables considered in the numerical experiments. 



o 

o 



0 

1 





^ o 




* 


""o D 

N 1 

+ +> 




o 








rH N 


o 




N <3 






w 


<\ N 




s: 


^k> 


+ 



N 

<3 



O 

fd 



o 



fd 

I 



o 



fd 

I 



I 

! 



in 

I 

o 



vD 

I 

o 



o 



B 

iH 

n| 



in 


in 


in 


in 


in 


>1 

•H 

C 


CN 




eg 




eg 


•H 


o 


o 


eg 


VD 


O 




• 


• 


• 


• 


• 


C 


o 


o 


o 


o 


eg 


H 




I 



B 



M 



in 

ro 



O 


in 


o 


in 


O 


in 


o 


in 


O 


in 


o 


in 


O 


in 


o 


in 


o 


in 


o 


o 


in 




o 


eg 


in 




o 


eg 


in 




o 


(N 


in 




o 


eg 


in 




o 


o 


o 


o 


iH 


iH 


iH 


iH 


eg 


eg 


eg 


eg 


ro 


ro 


ro 


ro 










in 


o 


o 


o 


o 


O 


o 


o 


O 


O 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


iH 



24 



TABLE II. Experimental tests that were carried out using the values in Table 
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Figure 1. Distribution of the mean wind velocity (U in 
m sec“l) with height (z in meters) utilizing series 1 for 
wave number (k) equal, respectively to; 1) initial mean 
velocity profile; 2) 0.475 m“l; 3) 0.375 m-1; 4) 0.275 m“^; 

5) 0.175 m“^; 6) 0.100 m“l; 7) 0.050 
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Figure 2. Distribution of mean wind velocity (U in m sec 
with height (z in meters) utilizing series 2 for wave number 



(k) equal/ respectively/ to: 1) initial mean velocity pro- 



0.475 m”l; 

6) 0.100 m“^; 7) 0.050 m“^, 



file; 2) 



3) 0.375 m 



- 1 . 



4) 0.275 m 



-1 



5) 0.175m 



-1 
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Figure 3. Distribution of the wave-caused stress (t in 
0.08 m^sec”^) with height for wave number (k) equal to 
0.350 m”^, and utilizing series 3. Variation in is 
equal, respectively, to: 1) 2.025 m; 2) 0.675 m; 

3) 0.225 m; 4) 0.075 m; 5) 0.025 m; 6) infinity. 
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Figure 4. Distribution of the wave-caused stress (T in 
0.09 m^sec”^) with height for wave number (k) equal to 
0.200 and utilizing series 3. Variation in is 

equal, respectively to: 1) 2.025 m; 2) 0.675 m; 3) 0.225m; 

4) 0.075 m; 5) 0.025 m; 6) infinity. 
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by Stauffer (1973) 



These two curves are very similar to 



the curves for 


^1 


= 2.025 m. 


changes rapidly 


in 


such a way 


z^ = 0.025 m is 


at 


least one 


than the case z 


1 


infinity . 


sensitivity of 


the 


solutions 



tion . 



As decreases, the curve 
that the maximum stress for 
order of magnitude greater 
These experiments show the 
to this lower boundary condi 



Series 4, 5, and 6, which use the finite difference 

estimate of the wind shear in the boundary condition, were 

carried out with three values of R . The solutions for 

m 

“5 “6 

R =10 and R =10 showed some significant differ- 
m m ^ 

-7 

ences while the solutions for R =10 were nearly identi- 

m ^ 



cal for those of R =10 

m 



-6 



Figures 5 and 6 show the mean 

-6 



wind profile for series 5 which employs R^ = 10 



-1 



These 

figures show that for k ^0.450 m * the wind profile is 
very close to the initial mean velocity profile^ Also 
note the large departures from the logarithmic profile at 
low levels for k = 0.1 m ^ or less. A disturbing feature 
of the solutions is the oscillation of the mean wind be- 
tween adjacent values of the wave numbers. 

Figure 7 shows the solutions for the mean wind for 
series seven. These solutions are the same as those of 
series 4 except that the grid size has been cut in half 
and the time step has been reduced to one-fourth. These 
solutions are generally similar to those in figures 5 and 
6, except that in the lower levels the wind is larger than 
the logarithmic value for k <0.42535 m These results 
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Figure 5. Distribution of the mean wind velocity (U in 
m sec~^) with height (z in meters) utilizing series 5 
fcr wave number (k) equal, respectively, to: 1) initial 

mean velocity profile; 2) 1.000 m“^; 3) 0.500 m“^; 

4) 0.450 m“^; 5) 0.400 m“^; 6) 0.350 m**!. 7) 0.300 m“^. 
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Figure 6* Distribution of the mean wind velocity (U in 
m sec“^) with height (z in meters) utilizing series 5 for 
wave number (k) equal, respectively, to: 1) initial mean 

velocity profile; 2) 0.250 m”!; 3) 0.200 m“^; 4) 0.150 m“^; 

5) 0.100 m“l; 6) 0.050 m“l. 
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Figure 7. Distribution of the mean wind velocity (U in 
m sec”^) with height (z in meters) utilizing series 7 
for wave number (k) equals respectively, to: 1) initial 

mean velocity profile; 2) 1.000 m“l; 3) 0.42535 m-1; 

4) 0.400 m-1; 5) 0.200 m"l; 6) 0.175 m"!; 7) 0.150 m"!. 
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Figure 8* Distribution of the mean wind velocity (U in 
m sec ) with height (z in meters) for wave number (k) 
equal to 1.000 m"^. Comparison is made where: 1) initial 
mean velocity profile; 2) series 4; 3) series 5 and 6; 

4) series 7. 
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Figure 9. Distribution of the mean wind velocity (U in 
m sec“^) with height (z in meters) for wave number (k) 
equal to 0.0400 m“^. Comparison is made where: 1) initial 

mean velocity profile; 2) series 4; 3) series 5 and 6; 

4) series 7. 
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Figure 10* Distribution of the mean wind velocity (U in 
m sec“^) with height (z in meters) for wave number (k) 
equal to 0.200 m*"^. Comparison is made where: 1) initial 

mean velocity profile; 2) series 4; 3) series 5 and 6; 

4) series 7. 
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Figure 11. Distribution of mean wind velocity (U in 
m sec with height (z in meters) for wave number (k) 
equal to 0.150 m . Comparison is made where: 1) initial 

mean velocity profile; 2) series 4; 3) series 5 and 6; 

4 ) series 7 . 
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Figure 12. Distribution of the mean wind velocity (U in 
m sec“^) with height (z in meters) utilizing series 8 for 
wave number (k) equal to 0.400 m“^. Variation in wave 
amplitude (a in meters) is equal, respectively , to : 

1) initial mean velocity profile; 2) 2.00 m; 3) 1.00 m; 

4) 0.50 m; 5) 0.25 m= .1/k; 6) 0.10m. 
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Figure 13. Distribution of mean wind velocity (U in 
m sec“^) with height (z in meters) utilizing series 8 for 
wave number (k) equal to 0.150 m“"^. Variation in wave 
amplitude (a in meters) is equal, respectively, to: 

1) initial mean velocity profile; 2) 2.00 m; 3) 1.00 m; 

4) 0.67 m = .1/k; 5) 0.50 m; 6) 0.10 m. 



40 



a 



(5.5) 



= 0.1 k 



-1 



In these figures, the amplitude ranges from 0.10 m to 
2.00 m. The smaller values of amplitude gives wind pro- 
files which are very close to the logarithmic profile. 

The larger wave amplitudes lead to increasingly larger 
departures from the logarithmic profile. This is to be 
expected since u is proportional to the square of the 
disturbance amplitude (equation 4.2). 

Figure 14 contains the disturbance potential tempera- 
ture as a function of z, from series 9 with three values 
of k. These curves show that the largest temperature 
amplitudes occur at about one meter in height; above these 
maxima the temperature decreases rapidly. The largest 
amplitude occurs for the smallest wave number. The 
numerical solutions showed that the stream function ampli- 
tudes were not significantly affected by the presence of 
temperature fluctuations. This was true for both T^ = 
+0.15 and T^ = -0.15. This behavior may be interpreted 
as follows. If it is assumed that the horizontal advec- 
tion is balanced by the vertical advection, then equation 
(2.17) becomes : 



0 



1 

U-c 




(5.6) 



This expression is used to estimate the last term in (2.19). 
It can be seen that the buoyance term will be much smaller 
than the vorticity advection term with the possible 
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Figure 14. Distribution of disturbance potential tempera- 
ture (0 in ®K) with height (z in meters) utilizing series 
9 for wave number (k) equal, respectively, to: 1) 0^275 m”^; 

2) 0ol25 m”l; 3) 0.050 m”^. 
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exception of the point where U = 
tions show that the presence of 
the conclusion that the buoyancy 



c. The numerical solu- 
such points does not change 
term is that unimportant. 
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VI 



CONCLUSIONS 



In this study the numerical results of Stauffer (1973) 
were extended with the addition of buoyancy effects and 
two different lower boundary conditions. The buoyancy 
effects have no appreciable affect on the stream function 
and wind stress fields; however^ the presence of buoyancy 
effects could affect the Reynolds wave stress through modi- 
fications in the mean wind profile. These effects were 
not considered in this study. 

In this thesis, emphasis was placed on modifications 
in the mean profile which arose from the Reynolds wave 
stress. The two new boundary conditions which were inves- 

3 u 

tigated involved two methods of estimating in one of 
the lower boundary conditions. This term was entirely 
neglected in the study by Stauffer (1973). In the first 
method the derivative was evaluated from the logarithmic 
wind profile evaluated at a false level (z^). This pro- 
cedure leads to very large wave stresses near the surface 
when became small. 

3 u 

In the second method, — was evaluated by finite dif- 

d z 

ferences. Many solutions were performed with this 
boundary condition involving different values of the wave 
number, convergence criteria, and the grid size. Those 
solutions showed that the wind profile was near the loga- 
rithmic profile for the larger wave numbers of order one. 
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Smaller wave numbers showed considerable fluctuation 
about the logarithmic profile in the lower layers with 
the tendency toward wind speeds larger than those given 
by the logarithmic profile. Figure 7, which employs the 
smallest grid size, shows this tendency clearly. However, 
further investigations with this profile suggest larger 
fluctuations if a finer convergence criteria is used. 

In general the computed wind profiles are consistent with 
observational data analyzed by Davidson (1974). 

Further investigations of this type must carefully 
consider the lower boundary condition which was shown to 
be very sensitive in this study. Perhaps a full two-dimen- 
sional treatment of the air flow over the wave surface is 
required. Also it is very important to obtain accurate 
expressions for the turbulent exchange coefficients at 
different locations near the wave surfaces. These studies 
should be coordinated with careful analysis of observa- 
tional data. 
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sensitive to the level where the velocity gradient is com- 
puted in the lower boundary condition, 4) disturbance 
potential temperature at lower levels increases inversely 
with wave number, and 5) the presence of temperature fluc- 
tuations has an insignificant effect on the Reynolds stress 
or on the stream function. The numerical solutions show 
general agreement with the observational analyses of wave 
modified wind profiles by Davidson (1974). However, these 
comparisons are not conclusive because of the sensitivity 
in the lower boundary condition. 
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